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Abstract. We present a brief review of Big Bang Nucleosynthesis (BBN). We dis- 
cuss theoretical and observational uncertainties in deuterium and helium-4 primordial 
abundances and their implications for the determination of important cosmological 
parameters. We present, moreover, some recent results on active-sterile neutrino oscil- 
lations in the early universe and on their effects on BBN. 

1 Introduction 

Big Bang Nucleosynthesis (BBN), as well known, is one of the solid pillars of 
the standard cosmological model. The theory predicts that relevant abundances 
of light elements, namely 2 H, 3 He, 4 He and 7 Li, have been produced during 
the first minutes of the evolution of the Universe. The predictions span about 
9 orders of magnitude and are in reasonable agreement with observations. The- 
oretical calculations are well defined and very precise. The largest uncertainty 
arises from the values of cross-sections of the relevant nuclear reactions. Theo- 
retical accuracy is at the level of 0.2% for 4 He, 5% for 2 H and 3 He and 15% for 
7 Li. However, comparison of theoretical results with observational data is not 
straightforward because the data are subject to poorly known evolutionary ef- 
fects and systematic errors. Still, even with these uncertainties, BBN permits to 
constraint important cosmological parameters and to eliminate many modifica- 
tions of the standard model, thus allowing to derive restrictions on the properties 
of elementary particles and, in particular, of neutrinos. 

In this paper, we briefly review the physics of BBN. In sect.|21we introduce 
the essential parameters and inputs. In sect.|2|we summarize the present situa- 
tion of observational data. In sect. 0] we discuss the determination of cosmological 
parameters. The last section is dedicated to BBN bounds on non-standard neu- 
trinos and, specifically, to BBN and neutrino oscillations 1 . 

2 The Physics of BBN 

To understand primordial nucleosynthesis, we must follow in detail the histories 
of nucleons in the early universe. This is usually done by using numerical codes 
(among which the Fortran code by Wagoner [2], updated by Kawano [3] has 

1 In this paper, due to space limitation, we will consider only selected topics. For a 
complete review of the BBN bounds on neutrinos see ref.Q 
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Fig. 1. Primordial light element abundances as predicted by standard BBN. The widths 
of the bands correspond to theoretical uncertainties. 



become a standard tool). However, the main features of Big-Bang nucleosynthesis 
can be obtained by simple analytic arguments described in this section (see 
|4I5I6I7I8I§| for details). 

Primordial nucleosynthesis occurs at temperatures T < 1 MeV, which are 
small with respect to nucleon masses. At these temperatures, the number of 
nucleons is simply equal to the initial baryon asymmetry of the universe. It is 
useful to describe this quantity in terms of the present baryon to photon ratio: 



V EE (N B - %)/iV 7 



(1) 



The parameter 77 is simply related to the baryon density of the universe, being 
f2 B h 2 = 3.7- 10 7 ry. 

The neutron to proton ratio is controlled by the weak processes 

n + e + <-> p + z/ e 
n + v e <-+ p + e~ 

n — > p + e~ + v e (2) 

which interconvert neutron and proton. When the temperature T of the universe 
is about 1 MeV, these reactions are fast enough to maintain neutron and proton 
in chemical equilibrium. The neutron abundance is thus given by: 



1 



l + exp(Zim/T + £ e ) 



(3) 



Big Bang Nucleosynthesis and neutrinos 3 



where Am = 1.29 MeV is the neutron-proton mass difference and £ e = Me/T is 
the dimensionlcss chemical potential of electron neutrinos (in standard BBN £ e 
is assumed to be negligible). 

When the temperature T drops below Tf = 0.6—0.7 MeV, the neutron-proton 
inter-conversion rate, 7w ~ G|T 5 , becomes smaller than the expansion rate the 
universe, H ~ \f g^G^T 2 , where counts the total number of relativistic degrees 
of freedom of the early universe. The deviation of g* from the standard value, 
,g* = 10.75, is usually described in terms of an equivalent number of massless 
neutrinos N v ^ 3 according to: 

0, = lO.75+J(JV„-3) . (4) 

For temperatures T < Tf , chemical equilibrium can no longer be maintained. 
The neutron abundance X n evolves only due to neutron decay, according to 
X n — X n (Tf) exp(— t/T n ), where r„ is the neutron lifetime. One should note 
that the "freeze-out" temperature Tf scales as Tf oc g\^ and thus is sensitive to 
the particle content of the early universe. The larger is g„ (or equivalently N v ), 
the earlier is the freeze-out of the neutron abundance, at an higher value, and 
hence, the larger is the 4 He abundance produced in BBN. 

When the temperature of the universe is equal to Tjv — 0.06 — 0.07 MeV 
neutrons and protons start to react each other to build up light nuclei. The exact 
value of Tv depends on the baryon to photon ratio r\. Only two body reaction 
are indeed important in BBN, such asp(n,7) 2 H, 2 H(p, 7) 3 He, 3 He(d,p) 4 He, etc. 
(see )10p. Deuterium must be produced in appreciable quantity before the other 
reactions can proceed at all. However, due to the large number of photons per 
baryon, photodissociation of deuterium is not suppressed until the temperature 
decreases well below the deuterium binding energy — 2.2MeV. Following 
jS], one can see that the temperature Tv below which deuterium production is 
favoured scales as T N ~ B d /(15 — In 7/). 

Once deuterium is formed, nucleosynthesis begins and light nuclei are pro- 
duced rapidly. Essentially all available nucleons are quickly bound into 4 He, 
which is the most tightly bound light nuclear species. In addition to He, sub- 
stantial amounts of 2 H, 3 He and 4 He are produced. No heavy elements (A > 8) 
are produced, due both to the fact that Coulomb-barrier suppression is very 
significant and to the absence of stable isotopes with A = 5 and A = 8. In fig. 1 
we show the light element abundances produced during BBN, as calculated by 
using the Kawano code 0, for ij between 10 -10 and 10 -9 . The calculation of 4 He 
abundance includes small corrections due to radiative processes at zero and fi- 
nite temperature, non-equilibrium neutrino heating during annihilation, and 
finite nucleon mass effects |11I12) . 

Theoretical predictions are affected by uncertainties at the level of 0.2% for 
4 He, 5% for 2 H and 3 He and 15% for 7 Li. These uncertainties are due to un- 
certainties in the weak rates (which are "normalized" to the measured neutron 
lifetime r n = 885. 7± 0.8 s, see |3 for details) and in the values of the relevant nu- 
clear reaction rates. They have been estimated by montecarlo or semi-analytical 
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methods |1(JI13| . Recently the nuclear data have been re-analyzed, leading to 
improved precision in the abundance predictions |14ll5ll5| . 

3 Observational Data 

The abundances of light elements synthesized in the Big Bang have been subse- 
quently modified through chemical evolution of the astrophysical environments 
where they are measured. The observational strategy is to identify sites which 
have undergone as little chemical processing as possible and rely on empirical 
methods to infer the primordial abundances. For example, measurements of deu- 
terium are made in in quasar absorption line systems (QAS) at high redshift; if 
there is a "ceiling" to the abundance in different QAS then it can be assumed 
to be the primordial value. The 4 He abundance is measured H II regions in blue 
compact galaxies (BCGs) which have undergone very little star formation. Its 
primordial value is inferred either by using the associated nitrogen or oxygen 
abundance to track the stellar production of helium, or by simply observing the 
most metal-poor objects. Closer to home, the observed uniform abundance of 
7 Li in the hottest and most metal poor Pop II stars in our Galaxy is believed to 
reflect its primordial value. (We do not consider 3 He whose post-BBN evolution 
is more complex.) 

As observational methods have become more sophisticated, the situation has 
become more, instead of less, complex. Relevant discrepancies, of a systematic 
nature, have emerged between different observers. In the following, we give a 
brief summary of the present situation for deuterium and helium-4 (looking from 
outside by a non-expert). We refer to 7 for a more complete and up-to-date 
discussion. 

Deuterium 

Post-BBN evolution of deuterium (D) is simple. Deuterium is burnt in stars 
producing 3 He. Any deuterium measurement provides thus a lower limit for the 
primordial D abundance and an upper limit for baryon density of the universe. 

In recent years, measurements of deuterium have been made in quasar ab- 
sorption line systems (QAS) at high redshift. These systems are presumably not 
contaminated by stellar processes and thus the observed deuterium should be 
close to the primordial one. Since deuterium isotope shift corresponds to velocity 
of only (—82) km/sec, clearly QAS with simple velocity structure are needed for 
reliable determinations. Moreover, ionization corrections, possible "interlooper" 
etc. further increase systematic uncertainties. 

In tab. n we give the results of recent deuterium determinations in QAS 
|17I18I19I2(JI21I22I23| . An estimate of primordial deuterium abundance can be 
obtained from the weighted mean of data in tab.^ R should be noted, however, 
that the dispersion among the different determinations is not consistent with 
errors in the single measurements (see [21] for detailed discussion). We will 
use, m the following, the value D/H = 2.78±g;f| x 10 given in [zj, which is 
the weighted mean of the log D/H values given by |17I18I19I2U<5T] . The quoted 
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error is the la error in the mean, given by the standard deviation of the five 
log D/H values divided by >/&. This error is used instead of the usual error in 
the weighted mean, in order to take into account the "anomalous" dispersion of 
deuterium data. 

Table 1. Deuterium abundance in quasar absorption line systems at high red-shift 
(see |21| for details). 



z 2.504 [TU 3.572 [H] 2.536 [E| 2.076 [201 2.526 |2] 3.025 a G2 

10 5 (D/H) 3.98t° lr 3.25 ± 0.3 2.54 ± 0.23 1.65 ±0.35 2A2±° ]f 5 3.75 ± 0.25 



a This system was first analyzed by with the result D/H = (2.24 ± 0.67) x 10~ 5 . 
The quoted value is from |22) . 



Helium-4 

As a result of stellar processing, 4 He is produced, increasing its abundance above 
the primordial value, together with "metals" , such as nitrogen, oxygen and other 
elements heavier than 4 He, which are not produced in BBN. The observed 4 He 
abundance provides thus an upper bound to the primordial one, Y p . 

Helium observations are done in H ff regions in blue compact galaxies (BCGs) 
which have undergone very little star formation (at present ~ 100 H II regions 
have been studied for their helium content). In order to infer the primordial 
value Y p , one extrapolates to zero metallicity [Z = 0) the observed relation 
between helium (Y) and metals (Z). This is usually done using by nitrogen 
(N) or oxygen (O) as metallicity tracers. Alternatively, one can simply average 
helium abundances in most metal poor objects. 

The present situation is that independent determinations of Y p have a statis- 
tical errors at the level of 1 — 2% but differ among each others by about ~ 5%. 
In particular, by using independent data sets, Olive and Steigman |24| and Olive 
et al. [251 have obtained Y p = 0.234 ± 0.003, while Izotov et al. and Izotov 
and Thuan [27| have found Y p = 0.244 ± 0.002. 

The discrepancy between different determinations is possibly related to differ- 
ent description of the complex physical processes acting in H II regions. Several 
sources of systematic uncertainties may, in fact, affect the helium determination 
at a level comparable or larger than the reported statistical errors, like e.g. the 
evaluation of the ionization correction factor (which is related to how much neu- 
tral helium is in the object under scrutiny), of the temperature correction factor, 
of underlying stellar absorption, etc. (see |28j for a review). 

As discussed in the next section, it is extremely important to have a better 
determination of the 4 He primordial abundances and a reliable estimate of the 
total (statistical + systematic) associated error. For our estimates, we will use 
the central values for Y p reported above (Yp = 0.234 and Y p = 0.244), and the 
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error estimate AY p = 0.005, which is obtained e.g. in |7I8| from the dispersion 
of the various Y p determinations. 

4 Cosmological parameters from BBN 

The deuterium abundance D/H = 2.78^0 38 x 10~ 5 can be used to determine the 
baryon density of the universe. As discussed in |21j . the quoted value corresponds 
to t) = 5.9 ±0.5 x 10~ 10 (in standard BBN) or, equivalently, to f> B h 2 = 0.0214 ± 
0.0020. The error budget is dominated by the observational uncertainties which 
are about a factor 3 larger than uncertainties in theoretical prediction. 

The obtained value for ft B h 2 has to be compared with independent deter- 
mination of the baryon density of the universe. In particular with the result 
ft-gli 2 = 0.0224 ± 0.0009 given in [23 which is obtained from a combined fit to 
the cosmic microwave background (CMB) and large scale structure (LSS) data. 
The agreement of these two independent determinations is extremely important 
because they rely on completely different physical phenomena (which occurred 
at different time during the evolution of the universe). We note that CMB (and 
LSS) are presently more accurate than BBN in determining the baryon density 
of the universe. 

The value of r\ deduced from deuterium can be used, in standard BBN, to 
predict the abundance of the other elements and to compare with observations. 
Following [21], one obtains Y p = 0.2476±0.0010, 3 He/H = 1.04±0.06 x 1(T 5 and 
7 Li/H = 4.5 ±0.9 x 10 -10 . It is evident that there is tension between the quoted 
values and the observational results. The "predicted" abundance for 4 He is higher 
than the "high" helium value of Izotov et al. |26I27| . Moreover, the "predicted" 
7 Li abundance is a factor 2-3 larger with respect to the present observational 
results jSO]. The origin of these differences has to be clarified. They could be 
due to systematic errors in the measurements or to evolutionary effects (e.g. 7 Li 
depletion) or they could be a real indication for non-standard effects in BBN. 

In particular, the present D and 4 He data seems to favour an equivalent num- 
ber of neutrino families N v < 3. In order to understand the present situation, 
it is useful to combine the deuterium value D/H = 2.78j^'gg x 10~ 5 , with the 
"low" helium abundance, Y p = 0.234 ± 0.005, or with the "high" helium abun- 
dance, Yp = 0.244 ±0.005. The error AY p = 0.005 is the "estimated" systematic 
error in 4 He measurements (see above). If we fit these data in the plane (rj,N v ) 
following J2], we obtain the bound N v = 2.3 ± 0.5 (la) in the first case, and 
N v — 2.8 ± 0.5 (la) in the second. In both cases, the central values are below 
three, even if the errors are large enough to allow for the standard value N u = 3. 

The described results clearly indicate that a large number of effective neutri- 
nos is disfavoured. One can conclude, in principle, that an upper bound on the 
number of extra neutrinos, SN V = N v — 3, is 8N V < 0.3. It is clear, however, that 
the situation is quite delicate. The error AN V — 0.5 is completely dominated 
by systematic error in 4 He measurements. For this reason, we believe that, at 
present stage, a more safe upper bound on the number of extra neutrinos is 
8N V < 1. Hopefully in the near future we will be able to derive a stronger limit. 
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Other physical parameters which can be bounded by BBN are the chemical 
potentials of different neutrino species, [i a where a — e,/i,T. The possible role 
of neutrino degeneracy in BBN was noted already in j^l] and then discussed in 
a number of papers. A non vanishing chemical potential for v e , or v T would 
increase the neutrino contribution to the energy density and can be described as 
an increase in N v . An additional (and dominant) effect exists for electron neutri- 
nos which directly participate to n-p interconversion reactions. A non vanishing 
(jL e would shift the equilibrium between neutrons and protons, see eq. ©, with 
large effects on light elements production. 

Several analysis have been made of the BBN limits on neutrino chemical 
potentials. A recent analysis which include also CMB data concludes: 



where £ a = fJ- a /T are the dimensionless chemical potentials. For further impli- 
cations and for a discussion of the case in which both N v and £ e are free to vary 
see 



5 BBN and neutrino oscillations 

Effects of neutrino oscillations on BBN are much different if only active neutrinos 
are mixed, if only one active and one sterile neutrino are mixed or if we consider 
the more "complete" case of mixing between three active and one sterile neutrino. 

5.1 Mixing between active neutrinos 

If neutrinos are in thermal equilibrium with vanishing chemical potentials, mix- 
ing between active neutrinos does not introduce any deviation from standard 
BBN results. The situation is more interesting if neutrinos are degenerate. In 
particular, it was shown recently that, for the mixing parameters which ex- 
plain the solar neutrino problem .'!2 (<5m 2 ol = 7.3 • 1CT 5 eV 2 and tan 2 so i = 0.4) 
and the atmospheric neutrino anomaly {Sm1 tmo = 2.5 • 10" 3 eV 2 and 
tan 2 # a tmo ~ 1), asymmetries in the muonic and/or tauonic neutrino sectors 
would produce, through oscillations, an asymmetry into the electronic neutrino 
sector. This means that, in presence of oscillations, the restrictive bounds on 
the chemical potential of electron neutrinos applies to all neutrino flavours. It is 
thus possible to obtain the restrictive bound: 



0.01 < £e < 0.2 

\Ur\ <2,6 



(5) 
(6) 



16*1 < 0.1 , 



(7) 



valid for any flavour |37|. 
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5.2 Mixing between one active and one sterile neutrino 

There are three possible effects on BBN created by mixing between active and 
sterile neutrinos. First is the production of additional neutrino species in the 
primeval plasma, leading to N v > 3. The second effect is a depletion of the 
number density of electronic neutrinos which results in a higher neutron freezing 
temperature. Both these effects lead to a larger neutron-to-proton ratio and to 
more abundant production of primordial deuterium and helium-4 (for the details 
see e.g. review £Q). If mixing between active neutrinos is absent the second effect 
would manifest itself only in the case of (y e — ^ s )-mixing, if we neglect relatively 
weak depopulation of v e through the annihilation D e v e — > v^.tV^^. 

The third effect is a generation of large lepton asymmetry due to oscillations 
between active and sterile species |3H]. However, this effect takes place only 
for very weak mixing, much smaller than the experimental bound and is not 
discussed in this paper. 

The problem of active-sterile neutrino oscillation is quite complex and has 
been discussed in many papers starting from 1990 (a large list of references can be 
found in ref. pQ). The problem was recently re-considered in |39| both analytically 
and by solution of the complete system of integro-differential kinetic equations. 
Earlier derived bounds have been re-analyzed and significantly different results 
have been found in the resonance case. 

The results of are shown in Fig. |2 The effect on BBN is expressed in 
term of variation of the effective number of neutrinos AN V . The upper panels 
are for the case of — v a (or v T — v s ) mixing, while the lower panels refer 
to the case of v e — v s mixing. The obtained results clearly depend on the sign 
of the mass differences. For positive 2 mass differences, 5m 2 > (left panels), 
sterile neutrino production occurs through non-resonant transitions. For Sm 2 < 
(right panels), one has instead resonant active-sterile transitions which result in 
much stronger bounds on the neutrino oscillation parameters. The solid lines in 
figEl correspond to numerical results, while the red dotted lines corresponds to 
analytic approximate results. It is evident that an observational bound on extra 
neutrinos much better than unity, say 5N V < 0.3, could give very restrictive 
limits on active-sterile neutrino mixing. Unfortunately, the present observational 
bound 8N V < 1.0 is not accurate enough to put relevant constraints. 

5.3 Three active and one sterile neutrinos 

It is practically established now that all active neutrinos are mixed with pa- 
rameters given by the Large Mixing Angle solution to solar neutrino problem 
(Sm 2 ol = 7.3 • 10~ 5 eV 2 and tan 2 6 so \ = 0.4) and by atmospheric neutrino data 
(8m 2 tmo — 2.5 ■ 10 -3 eV 2 and tan 2 # a tmo ~ 1)- Existence of fast transitions be- 
tween v e , z^, and v T may noticeably change BBN bounds on mixing with sterile 
neutrinos, expecially for small values of mass difference. In particular, due to 

2 In the notations of [39| . Sm 2 is positive if sterile neutrino is heavier than active 
neutrino, in the limit 9 — > 0. 
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Fig. 2. BBN bounds on active-sterile neutrino mixing. See f° r details 

oscillations between active neutrinos, a deficit of or v T would be efficiently 
transformed into a deficit of z/ e , leading to stronger bounds on active-sterile mix- 
ing. The effects of mixing between active neutrinos on the BBN bounds on a 
possible active-sterile admixture has been investigated in detail in |39|. 



6 Conclusion 

Comparison of BBN theoretical results with observational data is not straight- 
forward because the data are subject to poorly known evolutionary effects and 
systematic errors. Still, even with these uncertainties, BBN permits to constraint 
important cosmological parameters, like e.g. the baryon density il^h 2 , the effec- 
tive number neutrino families N„, the neutrino degeneracy parameters £ Q etc. 
The present bound on the number of extra neutrinos species SN U is about unity 
and is not accurate enough to put relevant constraints on active-sterile neutrino 
mixing. If this limit could be reduced in the next future, say to 8N V < 0.3, very 
restrictive limits on active-sterile admixture could be obtained. 
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